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Abstract 

We study a type I seesaw model of neutrino masses within the framework of A4 flavor symmetry. 
Incorporating the presence of both singlet and triplet flavons nnder ^4 symmetry, we construct 
the leptonic mass matrices involved in type I seesaw mechanism. We then constrnct the light 
nentrino mass matrix using the 3 cr valnes of neutrino oscillation parameters keeping the presently 
undetermined parameters namely, the lightest neutrino mass miightest, one Dirac CP phase 5 and 
two Majorana phases a, /3 as free parameters. Comparing the mass matrices derived using A4 
parameters as well as light neutrino parameters, we then evaluate all the A4 parameters in terms 
of light neutrino parameters. Assuming some specific vacuum alignments of A4 triplet flavon field, 
we then numerically evaluate all the free parameters in the light neutrino sector, using which we 
also find out the remaining A4 parameters. We then use the numerical values of these parameters 
to calculate baryon asymmetry through the mechanism of leptogenesis. We not only constrain the 
A4 vacuum alignments from the requirement of successful leptogenesis, but also constrain the free 
parameters in the light neutrino sector (miightest, a,/ 3 ) to certain range of values. These values 
can be tested in ongoing and future neutrino experiments providing a way to discriminate between 
different possible A4 vacuum alignments discussed in this work. 
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I. INTRODUCTION 


Understanding the dynamical origin of fermion masses and mixing has been one of the 
most widely investigated problems in particle physics during last few decades. Although 
the standard model (SM) of particle physics can explain the origin of mass through the 
Higgs mechanism, it can not provide a justihcation to the observed fermion mass hierarchy 
and mixing. Any attempt to gain an insight into such problems inevitably requires input 
from beyond standard model (BSM) physics. The observed patterns of quark masses and 
mixing still remains a puzzle and a significant number of research works have been done 
in order to understand their fundamental origin. The leptonic mass and mixing, after the 
discovery of tiny neutrino masses and their large mixing [ 1 ] have made the fermion mass and 
mixing problem even more puzzling. This is because the neutrino mass is found to lie at 
least twelve order of magnitude lower than the electroweak scale, and the pattern of leptonic 
mixing with large mixing angles is very different from quark mixing with small mixing 
angles. Recent neutrino experiments T2K [2], Double ChooZ [3], Daya-Bay [4] and RENO 
[5] have confirmed the earlier observations of tiny neutrino mass and large leptonic mixing 
and also measured the mixing parameters with more precision. Apart from more precise 
measurements of neutrino parameters, these experiments also provided strong evidence for 
a non-zero value of reactor mixing angle 6*13 which was thought to be (very close to) zero 
before. The best fit values of neutrino oscillation parameters that have appeared in the 
recent analysis of [ 6 ] and [7] are shown in table I and II respectively. 


Parameters 

Normal Hierarchy (NH) 

Inverted Hierarchy (IH) 

Am^i 

10-5eV^ 

7.50 

7.50 

1 

10-3eV^ 

2.457 

2.449 

sin^ 612 

0.304 

0.304 

sin^ 023 

0.452 

0.579 

sin^ 013 

0.0218 

0.0219 

Scp 

306° 

254° 


TABLE I: Global best fit values of neutrino oscillation parameters [6] 
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Parameters 

Normal Hierarchy (NH) 

Inverted Hierarchy (IH) 

Am^i 

10-5eV^ 

|Am§i| 

10-3eV^ 

7.60 

2.48 

7.60 

2.38 

sin^ 012 

0.323 

0.323 

sin^ 023 

0.567 

0.573 

sin^ 013 

0.0234 

0.0240 

Scp 

254° 

266° 


TABLE II: Global best fit values of neutrino oscillation parameters [7] 


Due to the absence of right handed neutrino, the SM can not explain neutrino mass 
through the conventional Higgs mechanism. Even if the right handed neutrinos are added 
by hand to the SM in order to allow a Dirac mass term, then the respective Yukawa couplings 
have to hne tuned to the level of 10“^^, which is highly unnatural. On the other hand, several 
BSM frameworks provide a natural origin of tiny neutrino mass by the conventional seesaw 
mechanism which can be broadly categorized into three types: type I [ 8 ], type II [9] and 
type III [10], all of which involve the introduction of additional heavy fermion or scalar 
particles into the SM. Similarly, several BSM frameworks have also been proposed in order 
to generate large leptonic mixing. Most of these frameworks introduce additional flavor 
symmetries, either discrete or continuous, into the SM. Prior to the discovery of non-zero 
6 ^ 13 , the neutrino oscillation data were in perfect agreement with some versions of p — r 
symmetric neutrino mass matrix. Out of four different neutrino mixing patterns that can 
originate from such a p —r symmetric neutrino mass matrix, the Tri-Bimaximal (TBM) [11] 
form of neutrino mixing received more attention in the literature. This particular mixing 
predicts the mixing angles as 623 = 45°, 612 = 35.3°, 6*13 = 0. However, since 6*13 = 0 has 
been ruled out by latest experimental data, one has to modify these /i — r symmetric or 
TBM type mass matrices in order to generate non-zero ^ 13 . Since the measured value of 
613 is much smaller compared to the other two mixing angles, one can still consider TBM 
mixing to be valid at leading order and can explain non-zero ^13 by incorporating small 
perturbations. This has led to several works including [12-18] within the framework of 
different BSM frameworks. 
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The TBM type mixing can be accommodated within several discrete flavor symmetry 
models [19]. Among them, the discrete group A 4 , group of even permutations of four objects, 
can reproduce the TBM mixing in the most economical way [20]. Necessary deviations from 
TBM mixing in order to generate non-zero 6^13 can also be explained within the A 4 model 
as shown by many groups including [13, 14]. Here we adopt the approach followed by the 
authors of [14] where they could accommodate all the neutrino parameters within a A 4 model 
by including additional flavons apart from those appearing in usual TBM analysis. Without 
assuming any specific realization of A 4 symmetry at leading order, we construct the most 
general neutrino mass matrix in terms of flavons assuming a type I seesaw framework. We 
then compare this mass matrix with the neutrino mass matrix constructed from the best fit 
values of neutrino parameters. To reduce the number of free parameters, we assume the A 4 
triplet flavons to acquire vacuum expectation values (vev) in a way which preserve Z 2 or Z 3 
subgroups of A4. This allows us not only to calculate the A4 flavon vev’s in terms of low 
energy neutrino parameters, but also to calculate the unknown neutrino parameters at low 
energy. These unknown low energy neutrino parameters include the lightest neutrino mass, 
which remains undetermined at neutrino experiments. This can be anywhere between 0 and 
the upper bound set by cosmological bound on the sum of absolute neutrino masses < 

0.23 eV [21]. The other unknown neutrino parameters are the leptonic CP violating phases: 
one Dirac CP phase and two Majorana CP phases. We determine these unknown neutrino 
parameters as well as the A 4 parameters for a particular A 4 vacuum alignment using the 
3cr values of experimentally measured neutrino parameters. This is the major advantage as 
well as motivation for choosing this particular approach as it allows us to determine all the 
neutrino parameters completely. 

After calculating the A 4 as well as unknown neutrino parameters for a particular A 4 vac¬ 
uum alignment, we then use the same set of parameters to calculate the baryon asymmetry 
of the Universe within the framework of leptogenesis. Apart from the origin of fermion mass 
hierarchy and mixing, the SM also fails to explain the observed baryon asymmetry. Although 
the BSM framework explaining the baryon asymmetry could be completely decoupled from 
the one explaining leptonic mass and mixing, it is more economical and predictive if the 
same model can account for both the observed phenomena. Leptogenesis is one such mech¬ 
anism which generates the observed baryon asymmetry by creating a leptonic asymmetry 
first and then converting it into baryon asymmetry through B + L violating electroweak 
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sphaleron transitions [22], According to the original proposal of Fukugita and Yanagida 
[23], this mechanism can satisfy all the Sakharov’s conditions [24] required to be fulfilled in 
order to produce a net baryon asymmetry. The out of equilibrium CP violating decay of 
heavy right handed neutrinos present in the type I seesaw mechanism can naturally produce 
the required lepton asymmetry. For a review of leptogenesis, one can refer to the review 
article [25]. Some interesting implementation of this idea within several BSM frameworks 
can be found in [26]. After calculating the baryon asymmetry in our model with type I 
seesaw and A 4 symmetry, we constrain the model parameters as well as different possible 
A 4 vacuum alignments by comparing our prediction with the observed baryon asymmetry 
given by [21] 

Yb = (8.58 ± 0.22) X 10"^^ (1) 

We hnd that certain combinations of A 4 vacuum alignments and light neutrino mass hierar¬ 
chy (normal or inverted) can not give rise to the observed baryon asymmetry. A few models 
can give rise to the observed baryon asymmetry for specihc values of the lightest neutrino 
mass and the leptonic CP phases. 

To allow experimental verification of the models, we also calculate, for each flavon align¬ 
ment, the effective neutrino mass ruee = (where U is the leptonic mixing matrix 

and rui are light neutrino masses) which can be probed in neutrinoless double beta decay 
experiments. We find that certain regions of parameter space can be ruled out by phase II 
of CERDA [27] experiment in future. 

This paper is organized as follows. In section II, we discuss our model of type I seesaw 
with A 4 flavor symmetry. In section III, we briefly discuss the mechanism of leptogenesis. In 
section IV we describe the numerical analysis adopted here and finally conclude in section 
V. 

II. TYPE I SEESAW WITH A4 FLAVOR SYMMETRY 

Type I seesaw [8] is the simplest possible realization of the seesaw mechanism where the 
SM held content is extended by three right handed neutrinos * — 1)2, 3) singlet under 
the SM gauge symmetry. This introduces additional terms in the Yukawa Lagrangian 

Cr = vt’e.Hi'i + i-'f C-'Af«4 + h.c. (2) 
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where £l = e)|^, H = (h°, h Y C is the charge conjugation operator. The resulting 

in 6 X 6 neutrino mass matrix after electroweak symmetry breaking is given by 


-M.= 


0 Md 

tT 


(3) 


Mjj Mjir 

where Mr, = yyV is the Dirac neutrino mass and v is the vev of the neutral component of 
SM higgs doublet. Assuming Mrr S> the light neutrino mass is given by the type I 
seesaw formula 


Ml = (4) 

Assuming the Dirac mass term to be at electroweak scale, one needs Mrr to be as heavy as 
10^“^ GeV in order to generate light neutrino mass of order 0.1 eV. The scale of right handed 
mass can be lowered by suitable fine tuning of the Dirac Yukawa couplings. 

We now briefly discuss the realization of type I seesaw mechanism that was presented 
in the work [14]. A 4 , the group of even permutations of four objects, is also the symmetry 
group of a tetrahedron. This group has 12 elements and four irreducible representations 
with dimensions rii such that — ^2. The characters of 4 representations are shown in 

table III. The complex number uj is the cube root of unity. These four representations are 
denoted by 1,1', 1" and 3 respectively. 


TABLE III: Character table of A 4 


Class 

Xi 

XV 

Xi" 

X3 

Cl 

1 

1 

1 

3 

C 2 

1 

OJ 


0 

G 3 

1 

a ;2 

u 

0 

G 4 

1 

1 

1 

-1 


Apart from the SM fields and three right handed neutrino fields required for type I 
seesaw, the model we are studying has five flavon fields, singlets under SM gauge symmetry, 
required to break the A 4 symmetry as well as to generate the desired structure of lepton 
mass matrices. The model also has an additional Z 2 symmetry to make sure the presence of 
only the desired terms in the Lagrangian. The transformations of all the fields in the model 
under A 4 x Z 2 symmetry are shown in table IV. The SU{2)l lepton doublets I = (G,^/ 4 ,G) 


6 






are assumed to transform as triplet 3 under A 4 whereas the SU (2 )l singlet charged leptons 
transform as 1 , 1 ", V respectively. The SU{2)l singlet right handed neutrinos 
transform as a triplet under A 4 . The SM higgs held H transforms like a singlet 1 under A 4 
whereas the SM singlet havon helds 4>Ei 4’Ni h; X) t/’ transform as 3, 3, 1,1', 1" respectively. 



H 

1 

l^R 

eR 

PH 

TR 

4>e 

<pN 

h 

X 

'P 

A 4 

1 

3 

3 

1 

1" 

T 

3 

3 

1 

T 

1" 

^2 

1 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

1 

1 


TABLE IV: Transformation of the fields under A 4 x Z 2 symmetry of the model 


After hxing the transformation of the model helds under the havor symmetry as well as 
the SM gauge symmetry, the Lagrangian for the lepton can be written as 


L D + XrTR){^^) + X^Hii/R + h.c. 

+^RRl^>R I 


(5) 


where A is the scale at which the havons acquire vev’s in order to break the A 4 symmetry. 
The dimensionless couplings CAr,c^,c^ and cp are, in general, complex. Due to the non¬ 
trivial Z 2 charge assignments, the havon held (pE couples only to the charged leptons. We 
can decompose the terms in the Lagrangian above into A 4 singlets using the A 4 product 
rules given in appendix A. Similar to the leptonic Lagrangian, one should also decompose 
the full scalar potential of the model into A 4 singlets and hnd out the vacuum alignments of 
the havon helds by minimizing the potential. Here we do not perform this detailed exercise 
which can be found elsewhere. Instead, we assume specihc vacuum alignments of the havon 
helds and study their phenomenological consequences. Assuming the triplet havon pE 
acquire vev as {Pe) = A.(l, 0, 0)^, we can write down the charged lepton mass matrix in the 
diagonal form as 


M, = 


A. 0 0 




0 A^ 0 


V 


( 6 ) 


0 0 A, 

where v is the vev of the SM higgs doublet as mentioned earlier. Similarly, decomposing the 


/ 
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term X^HiuR into A 4 singlets, the Dirac neutrino mass matrix can be written as 



Md= 0 0 Xnv 


( 7 ) 


^ 0 X]\fV 0 y 


The right handed neutrino mass matrix can receive contribution from four different flavons 
3^8 seen from the leptonic Lagrangian (5). Adopting the same notations as used 
by the authors of [14], we denote the triplet flavon vev’s as C 7 v( 07 v) = A{(f)a,(f>b,(f>c) and 
the singlet flavon vev’s as = Ar],c^{x) = = A'lp. Thus, in this notation 

0a, 6,0 h) Xy 0 are dimensionless and have absorbed the respective Yukawa couplings. In this 
notation, the right handed neutrino mass matrix can be written as 


/ 2 


\ 


|0a + h -|0c + 0 -|06 + X 


( 8 ) 


Mrr = Arr - 10 ^ + ^ |0b + X -\(t)a + 'n 

y -|0b + X -|0a + h |0c + 0 y 


We can use the Dirac and right handed Majorana mass matrices given by equations (7) and 
( 8 ) respectively in the type I seesaw formula given by equation (4) to construct the light 
neutrino mass matrix in terms of A 4 parameters. In the section on numerical analysis IV, 
we will discuss how these A 4 parameters can be computed in terms of the light neutrino 
parameters by comparing this mass matrix with the one constructed using light neutrino 
data. 

III. LEPTOGENESIS 

The observed Universe at present is baryon asymmetric, with the ratio of measured 
excess of baryons over anti-baryons to the entropy density is given by equation (1). If 
the Universe had started in a baryon symmetric way, three conditions must be satished in 
order to create a net baryon asymmetry. As pointed out hrst by Sakharov [24], these three 
conditions are (i) Baryon number violation, (ii) C and CP violation and (iii) Departure 
from equilibrium. Although the standard model satishes these conditions in an expanding 
Universe like ours, the amount of CP violation measured in the SM quark sector turns 
out to be too small to account for the entire baryon asymmetry of the Universe. This 
extra source of CP violation could be the leptonic sector which is not yet experimentally 






determined. Leptogenesis provides a minimal setup to connect lepton sector CP violation 
with the observed baryon asymmetry and also put indirect limits on these CP phases from 
the requirement of producing correct baryon asymmetry. 



FIG. 1: Right handed neutrino decay in type I seesaw models 


In a model with type I seesaw mechanism of neutrino masses, the CP violating out of 
equilibrium decay of the heavy right handed neutrinos into SM particles (as shown in figure 
1) can give rise to an asymmetry in the leptonic sector. This asymmetry is given by 

^ r(4 ^ L, + g-) - + 

< T{v% ^Li + H‘) + r(4 ^Li + H) ^ ’ 

If the three right handed neutrinos are assumed to be hierarchical such that M 2,3 ^ Mi, it 
is sufficient to consider the asymmetry produced by the decay of the lightest right handed 
neutrino z/^. Adopting the notations of [28], the lepton asymmetry resulting from the decay 
of can be written as 


e 


a 

1 


1 1 

87ru2 (mIMd)u 


J=2,3 


1 1 
8vru2 (mIMd)!! 




( 10 ) 


where u = 174 GeV is the vev of the Higgs doublet responsible for breaking the electroweak 
symmetry, 

g{x) = ^/x (1 + -(1 + x)ln^— 

\ 1 — a: X J 

and Xj = Mj/Mf. The second term in the expression for e" above vanishes when summed 
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over all the flavors a = e, fi, t. The sum over flavors is given by 


ei = 


1 1 


i=2,3 


The above leptonic asymmetry can be converted into baryon asymmetry as 


( 11 ) 


Yb = CK— (12) 

9* 

through electroweak sphaleron processes [22], Here, c is a measure of the fraction of lepton 
asymmetry being converted into baryon asymmetry and is approximately equal to —0.55. k 
is the dilution factor due to wash-out processes which erase the produced asymmetry and 
can be parametrized as [29] 


-K ~ V0.lKexp[-4:/{3{0.1Kf-^^)], for/sT > 10® 

r, for 10 < it' < 10® 


A"(lnA")0-® 

1 

2VK^ + 9 


, for 0 < a: < 10. 


(13) 


where K is given as 

^ _ Ti _ {MbMb)iiMi Mpi 

~ H{T = Ml) “ W 

Here Ti is the decay width of and H{T = Mi) is the Hubble constant at temperature 
T = Ml. Mpi is the Planck mass and the factor is the effective number of relativistic 
degrees of freedom at T = Mi which is approximately 110. 

The lepton asymmetry given by equation (11) is obtained by summing over all lepton 
flavors a = At very high temperatures (T > lO^^GeV) all charged lepton flavors are 

out of equilibrium and hence all of them behave similarly resulting in the one flavor regime 
of leptogenesis discussed above. However, at temperatures T < 10^^ GeV (T < lO^GeV), 
interactions involving tan (muon) Yukawa couplings enter equilibrium and flavor effects 
become important as discussed in details by the authors of [30]. Taking these flavor effects 
into account, the final baryon asymmetry can be written as 


y^2//ai?or _ 


-12 


[^29 


417 

5^’ 


+ ^l9 


390 


-m. 


Yr 


Sflavor 


B 


-12 


r e ^51 „ 

hr, 


I ^ 

+ eih 


/344 


-m. 




/344 


-m. 
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where €2 = e® + , m 2 = me + ifia = 

r]{nria) = 


Ml 


mn 


8.25 X 10-3eV 


-1 


. The function rj is given by 
-1 


0.2 X lO-^eV 


mn 


- 1.16 


For the calculation of baryon asymmetry, we first calculate the right handed neutrino 
mass spectrum by diagonalizing the right handed singlet neutrino mass matrix Mjiji as 


In this diagonal basis, according to the type I seesaw formula, the Dirac neutrino mass 
matrix also changes to 

mLR = m^Un (15) 

where is the Dirac neutrino mass matrix given by equation (7) in this model. 


IV. NUMERICAL ANALYSIS 


The light neutrino mass matrix can be constructed using the neutrino data of mixing 
angles and mass squared differences. The Pontecorvo-Maki-Nakagawa-Sakata (PMNS) lep- 
tonic mixing matrix is related to the diagonalizing matrices of neutrino and charged lepton 
mass matrices U^,Ui respectively, as 

f^PMNS = UjU,y (16) 


The PMNS mixing matrix can be parametrized as 


h^PMNS 


C 12 C 13 

-S 12 C 23 - ^25235136**^ 
S12S23 — Ci2C235l3e*'^ 


S 12 C 13 

C12C23 - Sl2S23Sl3e*'^ 
-C12S23 - Si2C23Sl3e*'^ 



S 23 C 13 


h^Maj 


C23C13 y 


( 17 ) 


where Cij = cos 6ij, Sij = sin6*jj and 5 is the Dirac CP phase. The diagonal matrix UMaj = 
diag(l, e*“, contains the Majorana CP phases a, f3. In the model we are working, the 

charged lepton mass matrix is diagonal given by equation (6) and hence Ui = /, the identity 
matrix. Therefore, in the diagonal charged lepton basis f/pMNS = U^. The light neutrino 
mass matrix can now be written as 


i-U = C/pMNsM^Wf/jMNS 


( 18 ) 
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where = diag(mi, m 2 , ma) is the diagonal form of light neutrino mass matrix. For 

the case of normal hierarchy, the three neutrino mass eigenvalues can be written as = 

diag(mi, y/m^ + Am^i, -^/mf + Am|]^), while for the case of inverted hierarchy (IH), it can 
be written as = diag('\/m| + Am^a — Am^;^, y/m| + Am^a, ma). If the light neutrino 

mass is given by the type I seesaw formula (4), then the right handed neutrino mass matrix 
can be obtained by inverting the type I seesaw formula 

Mrr = (19) 


Thus we have two expressions for Mrr-, one completely in terms of A/^ parameters given by 
equation ( 8 ) and the other in terms of light neutrino parameters and Dirac neutrino mass 
(19). The right handed Majorana mass matrix Mrr is complex symmetric and hence six 
independent complex parameters. Thus, comparing Mrr given by equations ( 8 ) and (19), 
we can write down six independent equations relating A 4 , parameters with the light neutrino 
ones. Solving these equations give us the six A^ flavon parameters 0a, 0b, 0c, X, "0 given 
in Appendix B. It should be noted that, in the earlier work [14], hermitian conjugate of the 
diagonalizing matrices were used instead of transpose in the definitions of light and heavy 
neutrino mass matrices given in equation (18) and equation (19). Since the Majorana mass 
matrices are complex symmetric instead of hermitian, using transpose of the diagonalizing 
matrices is more appropriate and hence we have adopted that convention. This change will 
also give rise to changes in the equations relating A 4 and neutrino parameters shown in 
Appendix B from the ones given in the earlier work [14]. This will also change the solutions 
of the equations relating different flavon vev’s for specihc A 4 triplet vacuum alignments, 
which we discuss below. 

The light neutrino mass matrix constructed from the neutrino oscillation parameters has 
four free parameters namely, the lightest neutrino mass mughtest, one Dirac CP phase 6 and 
two Majorana CP phases a, (3 which remain undetermined experimentally till now. If we fix 
these four parameters, then the six A 4 parameters can be determined using the equations 
given in Appendix B upto a factor F = We can reduce the number of these free 

parameters in the light neutrino sector, if we consider some specihc havon vev alignments 
like the ones discussed by the authors of [14]. These specihc alignments of the A 4 triplet 
vev’s (0a, 0b, 0c) can keep the Z 2 or Z 3 subgroup of A 4 unbroken in the havon space. Similar 
to the work [14], here also we do our calculations for both Z 2 , Z 3 preserving and ^ 2,^3 
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breaking triplet vev patterns. These triplet vev patterns are listed in the table V. Each 


Subgroup 

A 4 triplet vev alignment 

Z 2 

(1,0,0),(0,1,0),(0,0,1) 

Z 3 

(-1,1,1),(1,1,-1),(1,-1,1),(1,1,1) 

Breaking 

(0,l,l),(l,0,l),(l,l,0),(0,l,-l),(2,l,l),(l,l,2),(l,2,l),(l,-2,l),(l,l,-2),(-2,l,l) 


TABLE V: A 4 triplet vev alignments that preserve as well as break a Z 2 or Z 3 subgroup of A 4 

of these vev alignments give rise to two equations involving 0^, </)(,, 4>c- Since (j)a, 06, 4>c are 
complex, these two complex equations give rise to four real equations which can be solved 
numerically to determine the four light neutrino parameters miightest, <5,«,/d- For example, 
the vev alignment (1,0,0) gives 0^ = 1,06 = 0,0c = 0 where 0a, 6,c are given by equations 
(Bl), (B2), (B3) respectively in appendix B. After cancelling out the common factor F 
from the last two equations, we are left with two complex equations involving only the light 
neutrino parameters which can be found out by solving these equations numerically. In all 
the cases of vev alignments listed in table V, the common factor F in (j)a,b,c can be cancelled 
out. After Ending the four neutrino parameters numerically, we calculate the remaining three 
singlet A 4 flavon vev’s p, y, 0. The common factor F can also be calculated numerically for 
each of these cases. 

We use the values of A 4 parameters determined for a particular triplet vev alignment 
in the right handed Majorana neutrino mass matrix given by equation (8). If we use the 
analytical expression for the factor F = then the mass matrix can be calculated 
upto a numerical factor A^. It should be noted that although the parameter F gets fixed 
numerically for each possible triplet vev alignment, we still have the freedom to choose either 
A AT or A/j/j. We choose \n, the Dirac Yukawa coupling of neutrinos, in such a way that the 
lightest right handed neutrino mass falls in the appropriate flavor regime of leptogenesis. 
The chosen values of A^ are listed in table VI. Using these numerical values of the model pa¬ 
rameters, we then compute the baryon asymmetry for all the cases of A 4 vacuum alignments 
discussed above. 

It should be noted that in the earlier work [14], the most general set of unknown neutrino 
parameters (mughtest, <5, a, 0) were not calculated. Rather, the authors of [14] considered 
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specific choices of Majorana CP phases. Here we compute the most general set of such 
solutions in agreement with the A 4 triplet vacuum alignments. 



1 Flavor 

2 Flavor 

3 Flavor 

>2 

1 

1.5 X 10"3 

1.5 X 10“® 


TABLE VI: Values of for different flavor regime 
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FIG. 2: Normal Hierarchy with triplet flavon vev alignment (1,0,0) 
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FIG. 3: Normal Hierarchy with triplet flavon vev alignment (0,1,0) 
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FIG. 4: Normal Hierarchy with triplet flavon vev alignment (0, 0,1) 
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FIG. 5: Normal Hierarchy with triplet flavon vev alignment (—1,1,1) 
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FIG. 6: Normal Hierarchy with triplet flavon vev alignment (1,1, —1) 
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FIG. 7: Normal Hierarchy with triplet flavon vev alignment (1, —1,1) 
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FIG. 8: Normal Hierarchy with triplet flavon vev alignment (1,1,1) 
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FIG. 9: Inverted Hierarchy with triplet flavon vev alignment (1,0,0) 
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FIG. 10: Inverted Hierarchy with triplet flavon vev alignment (0,1,0) 
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FIG. 11: Inverted Hierarchy with triplet flavon vev alignment (0,0,1) 
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FIG. 12: Inverted Hierarchy with triplet flavon vev alignment (—1,1,1) 
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FIG. 13: Inverted Hierarchy with triplet flavon vev alignment (1,1, —1) 
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FIG. 14: Inverted Hierarchy with triplet flavon vev alignment (1, — 1,1) 

1 1 1 (IH) 



0.02 0.03 

m3 (eV) 




0 0.01 0.02 0.03 0.04 0.05 


m3 (eV) 


1 1 ^ 1 

r 

_1_1_1_ 

__j_: 


0 0.01 0.02 0.03 0.04 0.05 


m3 (eV) 


FIG. 15; Inverted Hierarchy with triplet flavon vev alignment (1,1,1) 
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/ 

X 
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/ 

X 

X 

X 

X 

(-2,1,1) 
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X 

X 

X 

X 


TABLE VII: Summary of Baryon Asymmetry Results 


V. RESULTS AND CONCLUSION 

We have studied a simple beyond standard model framework with type I seesaw mech¬ 
anism of light neutrino masses where the symmetry of the standard model is enhanced by 
an additional A 4 flavor symmetry. The structures of leptonic mass matrices are dictated by 
the underlying A 4 flavor symmetry. We choose the A 4 flavon helds in such a way that the 
charged lepton mass matrix remains diagonal and the Dirac neutrino mass matrix takes the 
form shown in equation (7). The right handed neutrino mass matrix shown in equation (8) 
is constructed by taking contribution from both singlet and triplet flavon fields under A4. 
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FIG. 16: Effective neutrino mass for triplet flavon alignments (1,0,0) and (0,1,0) 
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FIG. 17: Effective neutrino mass for triplet flavon alignments (0, 0,1) and (—1,1,1) 
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FIG. 18: Effective neutrino mass for triplet flavon alignments (1,1, —1) and (1, —1,1) 


1 1 1 (NH) 


0 1 1 (NH) 


> 

CD 


E 


0.009 

0.008 

0.007 

0.006 

0.005 

0.004 

0.003 

0.002 

0.001 

0 



0 0.002 0.004 0.006 0.008 0.01 


> 

o 


E 



mi (eV) 


mi (eV) 


1 1 1 (IH) 


0.02 

0.01 



mg (eV) 


0 1 1 (IH) 


0.055 

0.05 

0.045 

0.04 

0.035 

0.03 

0.025 

0.02 

0.015 



mg (eV) 


FIG. 19: Effective neutrino mass for triplet flavon alignments (1,1,1) and (0,1,1) 
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FIG. 20: Effective neutrino mass for triplet flavon alignments (1,0,1) and (1,1,0) 
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FIG. 21: Effective neutrino mass for triplet flavon alignments (1, 2,1) and (1,1, —2) 
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FIG. 22: Effective neutrino mass for triplet flavon alignments (0,1, —1) and (2,1,1) 


1 1 2 (NH) 


1 -2 1 (NH) 


> 

CD 


E 


0.0055 

0.005 

0.0045 

0.004 

0.0035 

0.003 

0.0025 

0.002 

0.0015 

0.001 



> 

o 


E 


0.004 

0.0035 

0.003 

0.0025 

0.002 

0.0015 

0.001 

0.0005 

0 



0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 


mi (eV) 


mi (eV) 


0 1 -1 (IH) 


1 -2 1 (IH) 


> 

CD 


E 


0.055 

0.05 

0.045 

0.04 

0.035 

0.03 

0.025 

0.02 

0.015 

0.01 



> 

<D 


E 


0.05 
0.045 
0.04 
0.035 
0.03 
0.025 
0.02 
0.015 

0 0.01 0.02 0.03 0.04 0.05 0.06 



mg (eV) 


mg (eV) 


FIG. 23: Effective neutrino mass for triplet flavon alignments (1,1, 2) and (1, —2,1) 
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1 Flavor 

2 Flavor 

(4'aj (pc) 

NH (mi(eV)) 

IH (m3(eV)) 

NH (mi(eV)) 

IH (m 3 (eV)) 

(1,0,0) 

0.0059 

0.001-0.018 

— 

— 

(0,1,0) 

— 

0.001-0.01 

— 

0.019,0.025 

(0,0,1) 

0.004,0.0053 -0.0061 

0.004 - 0.008 

— 

— 

(-1,1,1) 

0.0032,0.0032,0.0037 

0.00165,0.0026 

— 

0.0001 

(1,1,-1) 

0.0025 - 0.0037,0.0054 - 0.0062 

0.007,0.009 

— 

— 

(1,-1,1) 

0.0057 - 0.0064 

0.008 - 0.013 

— 

— 

(1,1,1) 

— 

0.0013,0.005 

— 

0.039 

(0,1,1) 

— 

0.004,0.009,0.014 

— 

0.02 

(1,0,1) 

0.00025,0.00215 

0.012 - 0.02 

— 

— 

(1,1,0) 

— 

0.0072,0.009 

— 

— 

(0,1,-1) 

— 

0.003,0.012 

0.012 - 0.017 

0.008 - 0.01 

(2,1,1) 

— 

0.0051,0.0121 

0.017,0.018 

— 

(1,1,2) 

0.02 

0.002 - 0.009 

— 

— 

(1,2,1) 

— 

0.007 - 0.013 

— 

0.025,0.045 

(1,-2,1) 

— 

0.008,0.033 

— 

0.03 - 0.055 

(l,l,-2) 

0.032 

0.0028 - 0.0083 

— 

— 

(-2,1,1) 

0.003 - 0.006 

0.0012 - 0.004 

— 

— 


TABLE VIII: Values of lightest neutrino masses giving rise to correct baryon asymmetry 


We also construct the right handed neutrino mass matrix by inverting the type I seesaw 
formula where the light neutrino mass matrix is evaluated in terms of the light neutrino 
parameters. Comparing these two mass matrices allow us to express the A 4 flavon vev’s in 
terms of the light neutrino parameters. Out of the light neutrino parameters, we use the 
3cr global £t values of two mass squared differences and three mixing angles, leaving four 
free parameters: the lightest neutrino mass, one Dirac CP phase and two Majorana CP 
phases. Instead of working in this four dimensional free parameter space, we evaluate them 
numerically by assuming the A 4 triplet flavon vev alignments to take some specihc form 
which can preserve a Z 2 or Z 3 subgroup of A 4 . We also consider alignments which break 
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these subgroups for illustrative purposes. Evaluating these free parameters in the light neu¬ 
trino sector also allows us to determine all the A 4 parameters numerically. We then use 
these parameters to calculate the baryon asymmetry through leptogenesis. The appropriate 
flavor regime of leptogenesis is decided by the strength of Dirac neutrino Yukawa coupling 
which is chosen by hand. For example, we show the parameter space of light neutrino pa¬ 
rameters which satisfy a specihc Z 2 or Z 3 preserving A 4 vacuum alignment in the hgures 
starting from 2 to 15. Each point shown in the first three panels of these hgures correspond 
to those values of light neutrino free parameters (mughtest, <5, a,/?) which give rise to one of 
the specihc A 4 vev alignments listed in table V. The fourth panel in these hgures show 
the baryon asymmetry in the two havor regime for the region of parameter space shown 
in the hrst three panels. We also calculate the baryon asymmetry in one and three havor 
regime, the summary of which can be found in table VII. It can be seen from table VII that 
none of the models give correct baryon asymmetry in the three havor regime. This puts a 
constraint on the Dirac neutrino Yukawa coupling whose strength decides the havor regime 
of leptogenesis in the models we are studying. Referring to the table VI, one can put a 
lower bound ~ 1.5 x 10“^ below which none of the models discussed in this work give 
rise to successful thermal leptogenesis. In the two havor regime where ~ 10“^ — 10“^, 
successful leptogenesis is possible with normal hierarchy only for the models with A 4 triplet 
vev alignment (0,1, — 1) and (2,1,1). For inverted hierarchy, seven models corresponding 
to vev alignment (0,1, 0), (—1,1,1), (1,1,1), (0,1,1), (0,1, —1), (1, 2,1), (1, —2,1) give rise to 
correct baryon asymmetry. In the one havor regime, where Dirac neutrino Yukawa coupling 
is of order unity, all the models with inverted hierarchy give successful leptogenesis whereas 
half of the models with normal hierarchy give rise to the desired lepton asymmetry, as shown 
in table VII. 

We also list the values of lightest neutrino masses in table VIII for which successful 
leptogenesis can be obtained for the models mentioned in table VII. It can be seen that the 
lightest neutrino mass can vary from 10“^ eV to 0.055 eV depending on the model. The 
corresponding values of CP phases can be obtained from the hgures showing their variation 
with lightest neutrino mass. Thus we can not only constrain the A 4 vacuum alignment from 
the requirement of producing correct baryon asymmetry through leptogenesis, but can also 
restrict the free parameters in the light neutrino sector to certain range of values. Even if 
the observed baryon asymmetry could have a diherent origin not related to the neutrino 


sector, these models can go through serious scrutiny at those experiments sensitive to CP 
violating phases and absolute neutrino mass scale. To have an independent experimental 
probe of the parameter space studied in this work, we also calculate the effective neutrino 
mass rUee to which neutrinoless double beta decay experiments are sensitive to. We plot 
the effective neutrino mass as a function of lightest neutrino mass for the parameter space 
allowed by A 4 vacuum alignment in the hgures starting from 16 to 24. One can see that 
phase II of GERDA experiments in future with sensitivity all the way down to 58-74 meV 
should be able to rule out some region of the parameter space discussed in this work. It 
is interesting to see from the hgures 2-24 as well the tables VII-VIII that within the same 
unbroken subgroup of A 4 namely, Z 2 or Z 3 , simple permutations of A 4 triplet vev’s can 
give rise to signihcant changes in the neutrino parameters (mughtest, <5,«,/d) as well as other 
pbservables like baryon asymmetry and effective neutrino mass. 

A possible extension of this work could be to show explicitly the realizations of A 4 vacuum 
alignment discussed here by minimizing the full scalar potential. Another important aspect 
not considered in the present work is to study the effect of renormalization group evolution 
(RGE) on the low energy parameters in the light neutrino sector. Here we have used the 
low energy best ht values of light neutrino parameters to solve the havon equations. Since 
the havons acquire vev at a high energy scale, one should consider the evolution of light 
neutrino parameters under RGE from low energy to the high energy scale. In the present 
work we have ignored these effects and have used the same low energy values of neutrino 
parameters to solve flavon equations as well as to calculate baryon asymmetry. We leave a 
detailed calculation incorporating the effects of RGE to future work. 
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FIG. 24: Effective neutrino mass for triplet flavon alignment (—2,1,1) 





Appendix A: A 4 product rules 


A 4 is a discrete non-abelian group of even permutations of four objects. It is also the 
symmetry group of a tetrahedron. This group has four irreducible representations: three 
one-dimensional and one three dimensional which are denoted by 1,1', 1" and 3 respectively. 
Their product rules are given as 

10 1 = 1 

1 ' 0 1 ' = 1 " 

1 ' 0 1 " = 1 
1 " 0 1 " = 1 ' 

3 0 3 = 1 0 1' 0 1" 0 3 ^ 0 3 ^ 

where a and s in the subscript corresponds to anti-symmetric and symmetric parts respec¬ 
tively. Denoting two triplets as (ai,bi,Ci) and {a 2 ,b 2 ,C 2 ) respectively, their direct product 
can be decomposed into the direct sum mentioned above as 

1 ^ Uia 2 -|- biC 2 -|- Cib 2 

1^ ^ C 1 C 2 -|- aib2 -|- bia2 
^ bib 2 -|- Cia 2 -|- aiC 2 

3s ^ (2aia2 — biC2 — Cib2, 2ciC2 — aib2 — bia2, 2bib2 — aiC2 — Cia2) 

3a ^ (biC2 - Cib2, aib2 - bia2, Cia2 - aiC2) 


Appendix B: A 4 Flavon vev’s 


0a — - (^ 12^13 ~ (e*^Ci2C23Sl3 “ S12S23) (c23Sl2 + e*'^Ci2Sl3S23)) 

mi ^ ^ 

P 

H-(C13S12 + (e*^C23Si2Si3 -|- C12S23) (C12C23 ~ e*^Si2Sl3S23)) (bl) 

m 2 

+ —(s?3 - e^*‘^Ci3C23S23) 
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= - ((e*^Ci2C23Sl3 - Si 2S23)^ + Ci2Ci3(c23Si2 + e*‘^Ci2Si3S23)) 

mi 
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-((e*'^C23Si2Sl3 + Ci2S23)^ + (C13S12(—C12C23 + 6**^512513523))) 

m 2 


(B2) 


0c — - ((e*^Ci25i3523 + 623532)^ + 612613(6*^612623513 — 5i2523)) 
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(2m3Ci25i2523(6^*"mi - m 2 ) (633 + 6 **^ 623513 )) 
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